Controlling Microscopic Friction through Mechanical Oscillations 
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We study in detail the recent suggestions by Tshiprut et al. [Phys. Rev. Lett. 95, 016101 (2005)] 
to tune tribological properties at the nanoscale by subjecting a substrate to periodic mechanical 
oscillations. We show that both in stick-slip and sliding regimes of motion friction can be tuned and 
reduced by controlling the frequency and amplitude of the imposed substrate lateral excitations. 
We demonstrate that the mechanisms of oscillation-induced reduction of friction are different for 
stick-slip and sliding dynamics. In the first regime the effect results from a giant enhancement 
of surface difi^usion, while in the second regime it is due to the interplay between washboard and 
oscillation frequencies that leads to the occurrence of parametric resonances. Moreover we show 
that for particular set of parameters it is possible to sustain the motion with the only oscillations. 

PACS numbers: 68.35.Af, 62.25.-l-g, 05.40.-a 46.55.-|-d 
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INTRODUCTION 

In the emerging field of nanoscale science, major ef- 
forts have been devoted to form a deep, predictive un- 
derstanding of the mechanisms involved in fundamental 
tribological processes, envisaging alternative innovative 
solutions to control friction [IJ . The ability to manipulate 
frictional forces is extremely important for many techno- 
logical applications. One may wish to reduce or enhance 
friction, eliminate the high-dissipative chaotic and stick- 
slip regimes of motion, and instead, to achieve smooth 
sliding. Early stages of motion and stopping processes in 
micromechanical devices and computer disk drives, which 
exhibit chaotic stick-slip, pose, for example, a real prob- 
lem. The difficulties in realizing an efficient control of 
friction are related to the complexity of the task, namely 
dealing with systems with many degrees of freedom under 
a strict size interfacial confinement. The achievements of 
nanotechnology have opened new perspectives and new 
basic scientific questions in this direction, where novel 
local probes give access to the study of friction at the 
atomic level. 

A novel approach for tuning frictional response, which 
has recently attracted considerable interest d-Q, is the 
mechanical control of a system, via externally imposed 
normal or lateral vibrations of small amplitude and en- 
ergy to the sliding interface. Manipulations by mechani- 
cal excitations, when applied at suitable frequencies and 
amplitudes, can drive particles, which are close to the 
interface, out of their potential energy minima, thus in- 
creasing considerably surface diffusion and mobility, and 
reducing friction. Through the dynamical stabilization of 
desirable regimes of motion (e.g., smooth sliding), which 
would be otherwise unstable, substrate mechanical os- 
cillations can reduce high dissipative sliding behaviors 
(stick-slip). 

Since the idea is not to change the physical proper- 



ties of mating interfaces, flexibility and accessibility are 
main relevant features of this approach. Controlling fric- 
tional forces has been traditionally approached by chem- 
ical means 0, [13] , usually by supplementing base lubri- 
cants with friction modifying additives. However, stan- 
dard lubrication techniques provide prescribed tribolog- 
ical properties which fit a certain range of parameters 
of the sliding objects (such as the load) and the ex- 
ternal forcing, and are expected to be less effective at 
the nanoscales. Here, on the contrary, frictional prop- 
erties can be tuned continuously by the frequency and 
the amplitude of the out-of-plane or in-plane vibrations 
of the sliding surfaces. The predicted effects should be 
amenable to atomic force microscopy (AFM) tests 111 Il2 1 
or applying ultrasound to the sample [isl [l3 | , or even in 
studies of contact mechanics of a probe interacting with 
oscillating quartz crystal microbalance surfaces 15l-ll7[. 

In this paper, we extend the recent results 18[ concern- 
ing the ability to modify tribological properties at the 
nanoscale by subjecting a substrate to lateral mechan- 
ical oscillations. The investigation is limited to a one- 
dimensional (ID) Tomlinson-like model. We show that 
microscopic friction can be suitably tuned and reduced 
by controlling the resonant frequencies and amplitudes of 
the imposed substrate lateral excitations. Depending on 
the model parameters, the numerical simulations reveal 
the details of the observed different regimes of motion 
(smooth sliding, stick-slip, inverted stick-slip) and of the 
occurring transitions among them. Moreover we show 
that for particular sets of parameters values the motion 
can be sustained just with the oscillations. 



MODEL AND NUMERICAL METHOD 

From a theoretical point of view, it has been frequently 
shown that, despite the simplified level of description, 
phenomenological models of friction have revealed the 
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ability of mimicking the main features of the tribological 
behavior observed both in nanoscale experiments or in 
more complex MD simulation frameworks [l9| . 

The Prandtl-Tomlinson model is, probably, the most 
widely used in interpretation of tribological experiments 
due to its simplicity and its ability of accounting for the 
main physical features of atomic-scale friction (see, e.g., 
I2I |20|). In this mathematical description, a particle 
of mass m is coupled, via a harmonic spring of stiffness 
if, to a moving stage which slides at constant velocity 
vg- Moreover, the particle experiences a periodic (sinu- 
soidal) coupling V{x) — U sm.{2'iTx/h) to the substrate, 
where x is the current position of the particle. Introduc- 
ing the viscous damping coefficient 7, which results from 
the energy dissipation to phonons or other excitations 
in the substrate, the equation of motion for the particle 
becomes 



mx -\- rwyx 



dV{x) 
dx 



K{x~VGt). 



(1) 



If the pulling spring is stiff enough (i.e., K is larger 
than the effective substrate spring constant Kq = 
max [V"{x)]), then there is always a unique equilibrium 
position Xcq ~ vat for the particle. Consequently, at 
small driving velocity, the friction will be linear in vg- 
Things become more interesting once Kq exceeds K. 
In this situation, the potential energy landscape expe- 
rienced by the particle is characterized, at certain in- 
stances of time, by more than one stable position. The 
time dependence of the combined spring and substrate 
potential reveals that mechanically stable positions dis- 
appear at certain instances in time due to the motion of 
the spring. The system is locked in one of the minima of 
the potential energy landscape until the increasing elastic 
stress of the spring allows overcoming the barrier. After 
that, the slider is accelerated. The potential energy of 
the elastic stress is converted into the kinetic energy of 
the slip event and, eventually, dissipated into heat via the 
damping 7-term. Then, the particle rapidly drops to the 
next nearest metastable minimum, where it locks again. 
As a consequence, for sufficiently small driving velocity, 
the dissipated energy per sliding distance is rather inde- 
pendent of Vg, in agreement with many experimental re- 
sults. Obviously, this simple mechanical description can 
provide only qualitative interpretations of the underlying 
tribological processes. The Tomlinson model is described 
in detail in many surveys, e.g., in Refs. with a 

long list of applications to concrete physical systems (in 
particular to those concerning tip-based devices). 

In this work we consider a Tomlinson-like model as that 
recently introduced in ref. ^18!]- In particular, in order 
to study the effect of lateral vibrations on ID diffusion 
and friction, the model describes the diffusion and driven 
dynamics of an AFM tip (a single mass driven by an 
external guide) subject to a substrate, which oscillates in 
the lateral direction with amplitude A and frequency /. 




FIG. 1: (Color online) Averaged friction force dependence on 
the oscillation frequency externally imposed to the substrate. 
The behavior is shown for three different values of the driving 
stage velocity vq: 1.6 (solid), 8.0 (dashed), 16.0 (dotted). 
The lateral force Fl = K{x — xq) is normalized to the static 
friction force Fq = max[V' (x)], and the oscillation frequency 

/ to the system characteristic frequency /o = 1 



The potential of such oscillating substrate becomes 



V{x,t) = C/sin 



27r 

— (2;-|-^sin(27r/i)) 




(2) 



The motion of the tip in lateral direction is governed 
by the equation 

mx = Fsub ix,t)~K {x - xg) + T {t)-m^{x~Vsub), (3) 

where Fsub is the force due to the substrate potential 
V{x,t), and xg = wg*; Vsub = -2tt f A siii{2iT ft) rep- 
resents the velocity at which the substrate oscillates. 
In particular, the last two terms take into account the 
thermal and dissipation effects in the framework of the 
Langevin dynamics 2l|. The random force T{t) and the 
viscous damping term 7 are related by the fluctuation- 
dissipation theorem: 



{T{t)r{t')) 2jmkBTS{t - t') 



(4) 



with ks and T denoting the Boltzmann constant and the 
temperature respectively. These two added terms can be 
thought of as representing degrees of freedom inherent in 
a real physical system which are not explicitly included 
in the adopted simplifled mathematical description. In 
particular, we assume that the damping term in Eq. (|3]) 
is determined by the energy dissipation to phonons and 
electron- hole excitations in the substrate. As a result, 
this term is proportional to the relative velocity of the 
driven tip and the substrate. Under these conditions, 
the substrate vibrations cause a time-periodic (ac) force 
acting on the tip, Fac = mA{2TTf)'^ sin(27r/f). This force 
presents the effect of inertia. 
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In the strong (weak) dissipative regime, when 7 is much 
larger (smaller) than the characteristic vibrational fre- 
quencies of the system, the motion is overdamped (un- 
derdamped). The numerical integration is carried out 
using the well-known velocity Verlet algorithm. For ev- 
ery value of the parameters and after reaching the steady 
state, the system characteristics of relevant physical in- 
terest are measured. Natural units used to character- 
ize the model for length, time and energy are A, /is, 
and meV. If not stated differently, m = 1, U = 100, 
b = I, K = 12.63, T = 0, and 7 = 32 define our de- 
fault parameter values. For parameters close to these 
default values, the system can be considered to behave 
as overdamped. Under these parameter ranges, the pre- 
dicted effects, as for surface diffusion resonance frequen- 
cies, should be amenable, e.g., to atomic force microscopy 
(AFM) tests using shear modulation mode, for which 
we typically have m « lO^^^ ^ 10"^'^ kg, U ~ 0.25 
eV, and b ~ 0.4 nm. For these values of system pa- 
rameters, the characteristic frequency turns out to be 
/q — [l/b)y^U/m, w 1 — 10 fis^^ and velocity unit is 



RESULTS AND DISCUSSION 

In order to investigate the effectiveness of lateral 
mechanical oscillations and their ability of improving 
nanoscale tribological properties, we have performed nu- 
merical simulations in different ranges of model param- 
eters where the "natural" system (without oscillations) 
may exhibit both high dissipative stick-slip dynamics and 
smooth sliding behaviors. 

At a fixed value of the amplitude A{= 1), Fig. [TJ shows 
the influence of the substrate oscillation frequency / on 
the time-averaged friction force < Fl >=< K{x — xg) >, 
for three distinct values of the driving velocity vg, corre- 
sponding respectively to stick-slip dynamics for vq = 1.6 
and to smooth sliding for vq = 8 and vg — 16. The 
frequency / is varied adiabatically and, for each of its 
value, Fl is averaged over a sufficiently long simulation 
time after reaching the dynamical steady state. Depend- 
ing on the external stage velocity, the results presented in 
Fig. [T] exhibit two types of resonance minima of friction 
which are determined by different physical mechanisms. 
As we will show shortly the low frequency resonance, 
which is most pronounced for the low driving velocity 
(e.g., Vg — 1-6), results from the giant enhancement of 
surface diffusion induced by the oscillations of the sub- 
strate, as recently discussed in details in ref. 
the contrary, for higher driving velocities (e.g., vg 
16.0), sharp resonance minima in the average lateral force 
< Fl > are caused by the interplay between the system 
washboard frequency VG/b and the oscillation frequency 
/• 

An understanding of the effects of lateral substrate 
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FIG. 2: (Color online) Diffusion D of the tip as a function 
of the normalized oscillation frequency f / fo, fo = \f^' 
finite temperature {ksT /U = 0.01), for three different values 
of the oscillation amplitude A. 



vibrations on the frictional properties of the system at 
low/moderate driving velocities, comes from the study 
of the diffusion coefficient for the undriven tip as a func- 
tion of the oscillation frequency [l8|. For low external 
frequencies /, the tip apex is able to follow the slow mo- 
tion of the substrate, performing just small vibrations 
around the potential minima. Under this condition, if 
the thermal energy is sufficiently smaller than the height 
of the potential barrier, the probability to escape from 
the potential well is exponentially small. For the optimal 
matching of the oscillation frequency / and the ampli- 
tude A, the diffusivity D may be considerably enhanced 
[2^ . In this case, the tip approaches the maximum of the 
substrate potential at the end of half cycle of the surface 
oscillations, where the driving force acting on the tip di- 
minishes. Then, even a weak thermal noise splits the 
ensemble of tips into two parts that relax to the neigh- 
boring minima of the surface potential, and the resonance 
enhancement of diffusion is observed. A further increase 
of the frequency leads again to localized dynamics; in 
contrast to the case of low frequencies here the tip over- 
comes the potential barriers and oscillates between neigh- 
boring minima of the surface potential. The resonance 
enhancement of surface diffusion leads to a corresponding 
reduction of the mean friction force at / ~ 3.32 which is 
most pronounced for low driving velocities (see Fig. [l]for 
VG — 1-6). We found that the curves of Fig. [2] calculated 
for three oscillation amplitudes A are not changed quali- 
tatively with increase of temperature, the diffusion reso- 
nance peaks are still visible for higher temperatures but 
greatly broaden and damped. The driving force acting 
on the tip is proportional to the amplitude of oscillations 
and, as a result, we observe that the resonance frequency 
decreases with A. 

By systematically scanning the phase space of the 
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FIG. 3: (Color online) During the first 2 ms the tip is driven 
at vg = 1-6 and then detached from the guide. The substrate 
oscillates uninterruptedly over the whole simulation with / = 
3.32 and A = 1.44. Clearly, the tip motion is completely 
carried by the substrate vibrations, with an average velocity 
highlighted by the dashed line. 



model parameters, it is possible to tmie the character- 
istics of the external oscillating excitations, in terms of 
suitable choices of amplitude and frequency, to achieve 
very favorable sliding conditions, for which friction al- 
most vanishes or even Ff, < on average. For a negative 
value of the lateral friction force, the driving guide, mov- 
ing at speed vq, does not pull the tip apex anymore, but 
it is instead pushed forward by the tip itself. An example 
of such a sliding situation is shown in Fig. [31 where, once 
the external guide moving at vg = 1-6 is detached from 
the tip, the direct tip motion is fully sustained by the 
only sinusoidal oscillations of the substrate, at the dif- 
fusion resonance / — 3.32. Due to highly nonlinear and 
strongly fluctuating dynamics, the determination of de- 
sirable set of model parameters achieving this "negative" 
friction regime is not always a trivial task. 

For higher driving velocities, as those shown in Fig. [1] 
(i.e., Vg = 8.0; 16.0), the mechanism of reduction of fric- 
tion is essentially different, stemming from the interplay 
between the system washboard frequency VG/b and the 
oscillation frequency /. This resonant behavior of the 
driven tip is clearly displayed by sudden and sharp drops 
of < Fl > at specific oscillation frequencies Jr, which 
are given by the formula fn — VG/Nb, with N integer. 

For an in-depth investigation of these features, we have 
analyzed the detailed dynamics during the transition 
across a resonance frequency — wg/2& for the case 
of an external guide velocity vg — 15 (see Fig. Due 
to a dynamical combination of fast sliding and substrate 
oscillations, for frequencies f ^ fn, the lateral force time 
series shows an asymmetric sawtooth shape resembling 
that of the usual stick-slip but with very tiny amplitude; 
as soon as / = fu the lateral force stops to oscillate, 
meaning that the guide-tip separation remains constant. 
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FIG. 4: (Color online) Fl/Fq time-dependence passing 
through a resonance frequency. At A = 1.44 and vg = 15, 
the oscillation frequency is increased from / = 7.49 (0 — 2000 
Ais) to / = 7.50 = vg/2 (2000 - 4000 /is), up to / = 7.51 
(4000—6000 /is). The lateral force oscillations are almost com- 
pletely washed out at the resonance frequency. The system 
shows, across the resonance, a transition from a stick-slip-like 
to an inverted stick-slip-like regime. 



Finally, when / > //? the synchronization is lost and 
Fl starts again oscillating with doubled frequency, but 
showing an inverted stick-slip-like behavior. The reso- 
nance marks in general the transition from a negative 
to a positive phase difference between the characteristic 
washboard frequency VG/b and the oscillation frequency 
/, which is characteristic for parametric resonances. This 
effect can lead to significant reduction of friction in the 
sliding regime of motion, as shown by sharp drops in av- 
erage and instantaneous lateral forces presented in Figs. 
[T] and m It should be noted that this mechanism of re- 
duction of friction in the sliding state is very different for 
that discussed above for the stick-slip regime. 

During the adiabatic increasing procedure of /, it has 
been noted that switching to the resonance fa at differ- 
ent times may result in distinct values of lateral friction 
force because of the different initial conditions. Due to 
commensuration between the washboard and the 

frequency / at the resonance //?,, this behavior shows 
in Fig. [5] the occurrence of only a discrete set of al- 
lowed values for < F^ >; the evenly separation between 
these states is determined by the dimensionless constant 
Kb/F^. 



CONCLUSIONS 

In summary, the model demonstrates the effectiveness 
of imposed lateral excitations for the reduction and con- 
trol of friction at the nanoscale under different regimes 
of motion (from stick-slip to smooth sliding) of a mat- 
ing contact. For an optimal matching of the oscilla- 
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FIG. 5: (Color online) Lateral force vs. time in approaching 
the resonance frequency / = 7.50 = vg/'2- Switching to the 
resonance at different times may result in a discrete set of 
distinct and evenly separated values of Fl due to different 
initial conditions. 

tion frequency and amplitude, this mechanical manip- 
ulation may drive, for example, an interfacial system out 
of its potential energy minima, thus increasing consid- 
erably surface diffusion and mobility, and reducing fric- 
tion. Unfavorable high dissipative stick-slip motions can 
be suppressed in favor of the dynamical stabilization of 
desirable smooth-sliding regimes. 

It has been shown that physical mechanisms of the 
oscillation induced reduction of friction are different for 
stick-slip and sliding regimes of motion. For stick-slip 
dynamics, at low driving speed, giant enhancement of 
surface diffusion takes place. In this case, the observed 
tribological response may turn to be quite sensitive to 
thermal effects, able to spoil the sharp resonant nature 
of this phenomenology. Differently, for smooth sliding 
regimes at high driving velocity, it is the interplay be- 
tween washboard and oscillation frequencies that deter- 
mines the occurrence of favorable parametric resonances. 
A very careful tuning of the model parameters has been 
proved to lead to the possibility of a directed motion sus- 
tained by the sinusoidal substrate oscillations only. 
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